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We found that the spectrum of the linear dispersion mode 
localized along the zero-mass line for variable mass 
Dirac fermions shows the nonzero spectral flow as a 
function of the external magnetic field B. This happens 
due to the Aharonov-Bohm effect, leading to the linear 
dependence of every eigenvalue on B. The nonzero spec- 
tral flow allows to use the magnetic field slowly vary- 
ing in time to control the energy of the linear dispersion 
mode. 



Illustration of the Aharonov-Bohm ring generated by the ZML. 
The particle is shown on the mass-square landscape 



1 Introduction Massless Dirac fermions are charge 
carriers in graphene |T| and topological insulators [2 |, two 
very popular subjects in contemporary physics. Geomet- 
ric and topological concepts (see, e.g., Refs. |[3]|4l) are of 
crucial importance for these systems. In particular, it was 
shown recently that nonzero spectral flow |5| is possible 
there |6,7|, which means a creation of electron-hole pairs 
from vacuum by adiabatically changing magnetic field |7|. 
The works |l6l|7l are rather formal. Here we give a simple 
physical example of a situation with nonzero spectral flow. 
It is based on our previous consideration |8| of electron 
states associated with zero-mass lines (ZML). These lines 
appear both in graphene on hBN and other substrates and in 
topological insulators (e.g., in (Hg,Cd)Te/HgTe/(Hg,Cd)Te 
quantum wells) |8 1. 

The peculiarities of dynamics of two-dimensional 
Dirac fermions localized in the vicinity of ZML were 
discussed in many different aspects, see fS\ and references 
therein. Classically it is described by the Hamiltonian 
function 

F(p,r) = Vp' + "^(i-)^ (1) 
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where r = {x, y} are the coordinates, p — {p^., py} are 
the corresponding momenta and m is the mass. The mass 
square landscape plays a role of the effective potential if we 
choose H{p,r)^ as an effective Hamiltonian. Thus every 
ZML gives rise to a one-dimensional channel, or waveg- 
uide, supporting the propagation of modes localized in its 
vicinity. Every mode is characterized by the transversal 
wave number related to the transversal energy. It turns out 
that there is a mode corresponding to zero transversal en- 
ergy. This quasi-one-dimensional mode has a linear disper- 
sion; i.e. its energy E of motion along the ZML is propor- 
tional to the momentum. 

Let us summarize the peculiarities of the linear disper- 
sion mode (LDM). Neglecting tunneling effects we can say 
that it is unidirectional, meaning that along a given ZML 
only the motion in a single direction is allowed. The Bohr- 
Sommerfeld quantization rule giving the semiclassical an- 
swer for the spectral series associated with the LDM reads 

(j) Prdr = 2^71 + (f)B + ■ ■ ■ , (2) 
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Figure 1 Illustration of the spectral flow associated with 
the magnetic field. The spectral flow is defined as the num- 
ber of crossings of the horizontal line E = taken with 
the sign of the derivative at the crossing point. In the figure 
spectral flow is equal to one. 



where the integration should be performed along the ZML, 
r is the length of this line counted from some point, pr 
is the longitudinal momentum, n is an integer, (t)B — 
for the curve without self-intersections is the Berry phase 
locally associated with the curvature and by ... we denoted 
some terms resulting from the displacement from the ZML 
(see 1 8 1). For the LDM we have pr = E, thus Q takes the 
form El ^ 2n{n + 1 /2) + . . ., where / is the length of the 
loop. 

It is easy to modify the quantization condition for the 
case of non-zero external magnetic field. Indeed, the closed 
ZML forms an Aharonov-Bohm ring in this case. Thus, an 
additional magnetic flux — BS through the loop bound- 
ing area S should be added on the right: 



El = 27r(7i + 1/2) +<P- 



(3) 



In this spectral series every eigenvalue linearly depends on 
the magnetic flux. If is changed (approximately) by a 
multiple of 27r, then the spectral series as a whole does 
not change at all, except that the numbering of eigenvalues 
is shifted. This shift in the numbering, which can equiv- 
alently be defined as the number of crossings of the lines 
E — E{<P) with the horizontal axis i? = 0, is called the 
spectral flow. Thus, the LDM gives a nontrivial physical 
example of nonzero spectral flow. Note that modes corre- 
sponding to nonzero transversal energy do not contribute to 
the spectral flow, because they always have a gap around 
zero. For the effective Hamiltonian corresponding to the 
LDM on a closed ZML, switching on a magnetic field 
whose magnetic flux is a multiple of 27r is equivalent to a 
gauge transformation on the ZML. Of course, this transfor- 
mation cannot be extended to the entire domain bounded 
by the ZML if the magnetic flux is nonzero. Thus, if we 
restrict ourselves to what happens in the vicinity of the 
ZML, we can adopt the point of view that the effective 
Hamiltonian does not change at all — we only choose a dif- 
ferent gauge. The above-mentioned gauge transformation 
can be used to define a one-dimensional vector bundle on 
the product of the ZML by the circle, and then the spec- 
tral flow coincides with the index of the elliptic operator 



that is obtained by adding the term d/dt to the effective 
Hamiltonian and acts on sections of this bundle, t being 
the coordinate on the circle. 

If the magnetic field slowly varies in time, the energy 
levels will still obey the dependence (j3]l. This permits con- 
trolling the energy of the LDM by means of a magnetic 
field: an increase in the magnetic field will lead to an in- 
crease in the energy and vice versa. 

Quantum-mechanically the dynamics of the Dirac 
fermion in the magnetic field is described by the Hamilto- 
nian 

H^(t{p-A) + a,m, (4) 

where a = {a^ , cTj,} are Pauli matrices and A is the vector 
potential. We consider a uniform magnetic field oriented 
along z-axis B = V x A. 

Let us consider a bent zero line 7 given by {x, y} — 
R(t), where r is a natural parameter, i.e. |R'(r)| = 1. In 
the vicinity of this line we introduce new variables r, <^ by 
the equality {x,y} = R(r) + ^n{T), where n is a unit 
normal vector on the curve at the point r. In the curvilin- 
ear coordinates r, ^ the Schrodinger equation H\I/ = E<F 
can be written as = E§, where if" = yjl - k^r)^^, 
kij) — — (R'(r)n'(T)) is the signed curvature, and 



H 



io-R'fr) d 1 , ,d 

■ — , - icrnir — 

Vl - efc(r) dT ^1 _ ^fc(r) ^'di 

ikcrn{T) 



a- A 



(5) 



2(l-^fc(r)) 

The Hamiltonian ([5]l can be simplified in the adiabatic 
approximation. We suppose that in the transverse direc- 
tion to 7 the wave function is localized at the scale much 
smaller than the radius of curvature; i.e. ^fc(r) <C L On 
the other hand, it is natural to assume that the curvature 
does not change significantly within this scale, which im- 
plies that ^fc'(r) /A:(r) <C 1. We will consider the magnetic 
field B corresponding to the magnetic length la = 1/ VB 
comparable to the radius of the curvature 1 /k{T). This pro- 
vides a magnetic flux of the order of unity through the area 
bounded by the closed curve 7. 

In the vicinity of the ZML we have A(r, ^) = 
A(t, 0) + A^(r, 0)^ + ... Inti-oducing the notation 
Ao(r) = A(r, 0) and using the equality Aq{t) = 
(AoR')R'(t) + (Aon)n(T) we find H ~ Hq + Hi, 



Ho = crR'(T)PV + crn{T)P^ + azTU, 



(6) 



ik 



Hi = (TR.\T)Ck{T)pr - -rTn(r) - aA'^ir, 0)C, (7) 

where = P^-(AoR'),Pt = -id/dT,P^ = ~id/d^~ 
(Aon). Let us choose a gauge A = {—By, 0, 0}. Then 
(tA'^{t,Q) ~ ~<JxBny{T). Replacing the operator p^- in 
(j7]i by the variable p^ we obtain symbols |9| Hq and Hi. 

In the adiabatic approximation the effective dynamics 
along the line 7 is described by the effective Hamiltonian 
Lb — Lq + Li. The symbols Lq and Li can be found from 
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general considerations ifTOl . As usual, Lq is an eigenvalue 
of Hq, i.e. Hqx — LqX- The symbol Li is given by 



dr dpr 



dpr dr 



(8) 



where (•)^ means the integration over ^. Introducing the 
notation = Pr ^ (AqR') and the function xo by the 
equality x — e''-^''"'^xo we reduce the calculation of Lq 
to the case without magnetic field [SJ. It gives Lq — —Pr- 
Similar to IH, we find from ^ 



Li = G(t) [k{T)Pr + B] 



fc(r) 



(9) 



where G{t) = 2(xoi'Cxo2)5 is the displacement from 
ZML, xoi = {nx - iny)xoi and xoi, X02 are real. The 
explicit expression for G(t) is given in Ref. |8|. Thus, the 

effective longitudinal equation LbiP = EiJ; reads 



Pr+G{T) k{T)Pr + B 



+ - ij^E^j, (10) 



Equation (lOi can be solved in the semiclassical approxi- 
mation. The solution is 



ip = exp —iEr + i ApfiR 



xexp(^i J G{T)[B-k{T)E]dT+^-j^ . (11) 

where 0^ = (f)/2is the Berry phase locally associated with 
the curvature |8|, dcj) = k{T)dT. Here (/> is none other than 
the angle between the tangent to the curve 7 and the x- 
axis, so that it gets an increment of 27r when going around 
7. Using the answer ( 11 1 it is easy to find the quantization 
condition for the case in which 7 is a closed curve: 



E I 



G{T)k{T)dT 



2tt {n + - 



B + j G{T)dT^ , (12) 



where the terms given by the circular integrals describe the 
small corrections to the length I of the contour 7 and the 
area S bounded by this contour owing to the displacement 
of the wave function from 7. Thus, the increase of the total 
flux by flux quantum corresponds to the replacement n — > 
n + 1. As a result, the spectral flow is equal to one, as 
illustrated in Fig. [T] 

Below we give a formal proof of this intuitively clear 
statement. Assume that the magnetic flux 



<P = b[s + j>G{T)dT 



is an integer multiple of 27r, = 27rm. Then it follows 
from the quantization condition ( [T2] i that the spectrum of 
the effective Hamiltonian Ls is the same as for B = 0. 
In other words. Lb and L = Lb\b=o are isospectral. In 
fact. Lb can be obtained from L by a gauge transforma- 
tion. Lb — U^^LU, where 



U = exp 



: J AodR -B J G{t) dr 



Consider the family of operators 

L{t) ^ {l-t)L + tLB, 



< t < 1. 



The spectral flow arising when the magnetic field B is 
switched on is just the spectral flow sf L{t) of this family. 
The latter can be expressed as the index of some elliptic 
operator. Namely, consider the first-order differential op- 
erator d/dt + Lit) on the product 7 x [0, 1]. We glue the 
endpoints of the interval [0, 1] together, so that the interval 
becomes the circle and this product becomes the torus 
7 X S^. The operator d/dt + L{t) becomes an operator 
with continuous coefficients on this product if we use the 
transformation U to glue together the values of the func- 
tions on7x{0} 7x{l}. One can readily verify that this 
operator on 7 x S"^ is elliptic, and by Q 

sfL(t) = ind(^^+£(t)^ , 

where 'ind' denotes the index of the operator Thus, the 
spectral flow is equal to the number of zero-modes of the 
operator d/dt + L(t) minus that of the adjoint operator, 
that is, one. 

To conclude, we have shown that a closed zero-mass 
line naturally generates the Aharonov-Bohm ring for the 
Dirac fermions if the external magnetic field is applied. 
For the linear dispersion mode the Aharonov-Bohm effect 
shifts the entire (sub)spectrum, leading to nonzero spec- 
tral flow. This effect can be used to control the energy of 
the linear dispersion mode by means of the magnetic field 
slowly changing in time. The spectral flow can be com- 
puted in topological terms as the index of an elliptic oper- 
ator on the two-dimensional torus. 
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